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Abstract 

We present the exact diagonalization of the Schrodinger operator corresponding to 
a periodic potential with N deltas of different couplings, for arbitrary N. This 
basic structure can repeat itself an infinite number of times. Calculations of band 
structure can be performed with a high degree of accuracy for an infinite chain and 
of the correspondent eigenlevels in the case of a random chain. The main physical 
motivation is to modelate quantum wire band structure and the calculation of the 
associated density of states. These quantities show the fundamental properties 
we expect for periodic structures although for low energy the band gaps follow 
unpredictable patterns. In the case of random chains we find Anderson localization; 
we analize also the role of the eigenstates in the localization patterns and find clear 
signals of fractality in the conductance. In spite of the simplicity of the model many 
of the salient features expected in a quantum wire are well reproduced. 
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Quantum wires represent the dreamed idea of make a conductor wire as small as a 
molecule. The idea that macroscopic devices we have been using for years can actually 
be built in nature at the nanoscale size goes back to Feynman and is the basis of all work 
currently carried out in the field of Quantum Electronics. This area of research is not only 
interesting in itself from the fundamental point of view but has also profound implications 
in applied physics and material sciences as hundreds of experiments are being carried out 
nowdays with a high degree of success. The purpose of this paper is to show that some 
of the ideas underlying the actual development of different types of quantum wires can 
actually be modelled in quite a simple manner using elementary quantum mechanics. 
This is specially important from our point of view as it provides a bridge between current 
fundamental research and basic concepts of quantum mechanics which are usually the 
subject of graduate standard programs. 

The aim to modelate a simple one-dimensional solid in order to study its band structure 
goes back to Kronig and Penney in the thirties [TJ but since then much work has been 
done along the lines of this first seminal reference. A sample of the variety of models that 
can be constructed with the same ideas can be found in [T2J . As the main interest of almost 
all of these authors was band theory, the techniques used in all these papers were mainly 



addressed to semiconductor physics |L4| . Much more recently a revival of the same models 
and techniques |TjJ has arisen as a consequence of the interest in truly theoretical and 
experimental one-dimensional physical systems from which quantum wires are just only 



one example Right now the research in molecular conductances, one-dimensional 



metallic rings and other devices of the same sort lends support to the idea of generalizing 
old methods yielding exact analitic solutions coupled to the use of desk-top computer 
algebra. 

Here firstly we shall solve analitically the band structure of an infinite periodic 
chain of delta potentials each one with a different coupling inside the primi- 
tive cell paying mostly attention to the mathematical aspects. After studing the band 
structure of the chain one can introduce randomness boosted by quantum fluctuations in 
order to account for localization. Fortunately the model seems rich enough to yield more 
information such as the generalization of the Saxon-Hutner conjecture W71 and even scal- 
ing exhibited as fractal behaviour of the conductance |I8| . The Paper will be organized 



as follows. In Section I we shall present the analytic solution of the periodic case of 
infinite different delta potentials. The band structure of this one dimensional Periodic 
Potential will also be briefly discussed. In Section II we turn our attention to the case 
of random arrays of one dimensional delta potentials with different couplings. Here we 
discuss the density of states using the functional equation method and classify the dif- 
ferent types of localization appearing when disorder is present. As we have been able to 
increase our understanding of localization to the extent of detecting universality effects, 
we shall entirely devote the Section III to the discussion of fractality in the conductance 
and the different checkings we have managed to perform in order to ascertain ourselves 
and try to convince the reader that the effect is present even in this extremely simple 
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model. We close with a Section of Conclusions. 



1 Periodic Array 



Let us consider an electron in a periodic one dimensional chain of atoms modelled by the 
potential constituted by an array of N delta functions each one with its own coupling ef , 
(i = 1, 2, ...N). After finishing the iV-array, the structure repeats itself an infinite number 
of times. The number of species N, can be arbitrarily large but finite. The case N=l is 
an old textbook exercise but may be convenient to be revisited |T!| for taking a full profit 
of our general results. The generalization can thus be followed in a more straightforward 
manner. The relevant primitive cell for N = 2 can be represented for the following set of 
wavef unctions: 



primitive cell 



= A x e ikx + B x e~ ikx < x < a 

yjf 2 (x) = A 2 e ik ^-^ + B 2 e- ik( - x -^ a < x < 2a 

* 3 (x) = e 2iQa [A ie ik{x - 2a) + B x e- ik{x - 2a) ] 



The matrix relating the amplitudes of the above wave functions for this N = 2 case can 
be written as: 



Aka 



( e" 
—ike tka 



ike~ lka 



—J- —T 



-ika 



^ika 

-ike tka 
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ike~ lka 



J 



It is trivial to generalize these two steps to the case of three species (i.e. N — 3). One 
can equally write the correspondent matrix in the form: 

/ 





g— ika 
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—ike lka 


ike~ lka : 
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(ik - ^) 

V a2 / 











gifca 


g— jfca 
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—ike lka 


ike~ tka j 



(2) 



where in all of the above cases, we have used the notation: 



= exp{iNQa} and Q G 



7T 7T \ 

Na'Na) 



(3) 
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and we shall also be using the length of each species, defined as: 

n 2 



n 9 

met 



(4) 



One can now proceed to the generalization of the matrix form for general number N of 
species just by defining the following 2x2 matrices: 

gifca g— ika ^ j J | 



E -[_ ike ika lke -ika) 5 Aj-f^_^ _^ + A)J- (5) 



The matrices for N = 2 and AT = 3 species given by (|T|) and @ can now be put in a 
more compact form with the help of the E and Aj as: 



E 



Ax 
E 



4x4 




Ax 


02x2 


E 


A 2 


02x2 


E 



(6) 



6x6 



It is now relatively simple to guess that the general form of a matrix for iV species 
must be written as: 



/ E 

02x2 



Ax 
E 







2x2 



02x2 

A 2 
E 

02x2 



02x2 

A3 02x2 



02x2 
02x2 

02x2 



\ 



02x2 

V^A N 







2x2 







2x2 



E 

02x2 



02x2 

A N _i 
E / 



(7) 



2Nx2N 



So far nothing very exciting has happened except that one can write the matrices in a 
compact, logic and generalizable way. And in fact without further steps the progress would 
have not been certainly remarkable. The real breakthrough arises when one realizes that 
one has to deal with the determinants equated to zero of these matrices in order 
to learn something about the band condition of this one dimensional N-species 
quantum periodic structure. Let us define the following function (e = ka) []: 



hA 



cos(e) + 



a \ sin(e) 



1 Notice that for negative energies, e takes pure imaginary values that we represent in the figures in 
the negative part of the spectrum 
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A quite simple computer algebra calculation shows that the determinant equated to zero 
of (HD, which belongs to the N = 2 case, can be written in terms of these functions as 

cos(2Qa) = 2h x h 2 - 1 (9) 

And the determinant of the N = 3 matrix given by (Q) can also be calculated to yield: 

cos(3Qa) = 4/ii/i 2 /i 3 - (hi + h 2 + h 3 ) (10) 

Below, the cases N = 4, 5, 6 and 7 are explicitely listed, using the generalized matrix (^) 
for each case and calculating the determinant equated to zero with the help of the h(e) 
functions (||). The result is: 

cos(4Qa) = 8 h x h 2 h 3 h 4 - 2 (h 4 h 2 + h 4 h 4 + h 2 h 3 + h 3 h 4 ) + 1 (11) 

cos(5Qa) = 16 h\h 2 h 3 h 4 h^ — 4 (h x h 2 h 3 + h x h 2 h 5 + h x h 4 h 5 + 

h 2 h 3 h 4 + h 3 h 4 h 5 ) + (hi + h 2 + h 3 + h 4 + h 5 ) (12) 

cos(6Qa) = 32 h\h 2 h 3 h 4 h^h^ — 8 (h x h 2 h 3 h 4 + h x h 2 h 3 hQ + hih 2 h^h Q + 

+ hihji b h§ + h 2 h 3 h 4 h 5 + h 3 h 4 h 5 h 6 ) + 2 (hih 2 + h x h 4 + h x h e + 

+ /12/13 + h 2 h 5 + /i 3 /i 4 + /i 3 /i 6 + h 4 h 5 + h 5 h e ) - 1 (13) 

cos(7<5a) = 64 h\h 2 h 3 h 4 h^h^h 7 — 16 [h\h 2 h 3 h 4 h^ + h x h 2 h 3 h 4 h 7 + 

+ h x h 2 h 3 h G h 7 + h\h 2 h$h§h 7 + hih 4 h b h & h 7 + h 2 h 3 h 4 h 5 h 6 + h 3 h 4 h 5 h 6 h 7 ) + 
+ 4 (hih 2 h 3 + h x h 2 h 5 + h x h 2 h 7 + h x h 4 h 5 + h x h 4 h 7 + h\h§h 7 + h 2 h 3 h 4 + 
+ h 2 h 3 h 6 + h 2 h 5 h 6 + h 3 h 4 h 5 + h 3 h 4 h 7 + h 3 h 6 h 7 + h 4 h 5 h 6 + h 5 h 6 h 7 ) - 

- (Tii + /i 2 + /is + /i4 + /is + /ie + /i?) • (14) 

We have been able to proof by induction that the general form of the iV species case 
can be given as: 

cos(iVQa) = S(e;ai,... , ajv) (15) 

• iV even 



i3(e;a!,... , ajv ) =2 Af - 1 ^/ ij ...(iV).../ ifc - 2*" 3 2)...fc fc + 

p p 

+ 2 V •^/,,,..(.V - 1)...//, - ...( 1)- 1 2^/;,..f2)...//, + (-l)f (16) 
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N odd 



B(e; a u ...,a N ) = ^ /,,...( \)...h h - 2 N ~ 3 ^ ^...{N - 2). 

p p 


..h k + 




+ 2 N ~ 5 £/,,,..(. Y -A)...h k - ... (-1)^ 2 2 £ ^...(3). 

P P 


..h k + 


(17) 


N—l 

(-1) — (ht + h 2 + h 3 + . . 


+ h N ) 





All what remains is to define the symbol hi...(M)...hk which means the sum of all 
possible products of M different hi 's with the following rule for each product: the indices 
must follow an increasing order and to an odd index must always follow an even index 
and reciprocally. 

The band structure provided by ([16]) and fllTD is not just exact but also extremely 
useful from the point of view of computer algebra calculations. In fact we have carried 
out various profiles for the curves provided for these conditions until iV=30 or more using 
just few seconds of a lap-top regular computer. The reason for that lies mainly in the 
systematic use of the form, products and combinations of the /i(e)-function defined by 
(H). As examples of what has just been said we list in Fig. [I] and || a series of band 

curves for large number of species and various values of the parameter f , defining the 

characteristic value of the /i(e)-function. One can observe also the unpredictable set of 
allowed bands which appear at low energies. This pattern increases its unpredictibility 
with the number of species. 

Once the band condition is known, one can write the distribution of electronic states 
in a very simple form. In one dimension the density of states per unit length of the chain 
for the nth band comes from, 



9n{e) 



2tt 4^ 



de(Q) 



dQ 



where the sum is over all the first Brillouin zone (1BZ) points Q with the same energy e. 
Due to the parity of cos(NQa) the number of points with the same value of e in the 1BZ 
is always 2, and provided that overlapping of neighbouring bands is not possible in this 
system, we can write the density of states as 



dQ{e) 



de 



(19) 



inside the permitted bands. From ([15]) a trivial calculation leads to, 

dB(e) 



G(e) 



g[e) ■ a 



[1 ~ B 2 ] 

Ntt 



de 



(20) 
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Figure 1: Band condition (detail in inset) and band structure (1BZ) for a periodic chain 
with six species in the primitive cell (~): 2, 2.5, 1, 1.5, -1, -2. 

where G(e) is the density of electronic states per atom. Fig. ^ shows some examples of the 
characteristic form of the distribution of states for different configurations of the primitive 
cell. 

2 Random Chains 

The structures one can observe in Nature hardly show a perfect periodicity. Even in 
the laboratory it is a difficult task to grow a crystal free of impurities, vacancies or 
dislocations. We shall now treat the presence of substitutional disorder in one dimensional 
delta-potential chains, that is we consider a chain of equally spaced deltas in which the 
sequence of different species does not obey a periodic pattern. This model has been mainly 
studied regarding the vibrational spectrum 0, paying less attention to its electronic 
density of states |3). For the purpose of studying quantum wires the relevant behaviour 
we want to analyze lies more in the latter than in the former physical property. 
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Figure 2: Band condition (detail in inset) and band structure (1BZ) for a periodic chain 
with ten species in the primitive cell (— ): 1, 2, 3, -1, -2, -3, 0.5, -0.5, 2, I. 




(a) Two species (f-): -0.3,-0.9 (b) Four species (f-): -3, -5, -6, -4 

Figure 3: Density of electronic states for different configurations of the primitive cell. 
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2.1 Energy gaps 



The Saxon and Hutner conjecture ]T?| was proved by Luttinger for the case of binary 



chains 0. We have been able to extend it to the general case. A detailed calculation 
following the line of Schmidt can be found in Appendix The result can be easily 
summarized as follows: the forbidden bands that coincide in different one species 
delta chains with couplings ef, . . . ,e 2 N are also forbidden levels in any infinite 
chain made up of deltas of the N types. This conclusion can be applied to both 
disordered and periodic chains. As can be seen in the calculations the result requires the 
interatomic distances of the chains involved to be constant and the same for all of them. 



2.2 Eigenenergies for finite disorder chains 

Let us calculate the allowed energy levels of a finite non-periodic chain of deltas with fixed 
end-points boundary conditions. The procedure to follow is the same as for the periodic 
case but imposing the vanishing of the wave function at the end-points of the chain which 
we locate at one atomic distance to the left of the first delta and to the right of the last 
one. The connection of the wave function throughout the different sectors leads us to 
a condition for the permitted energy levels in the form of the determinant of a matrix 
equated to zero. Again these matrices can be written in a generalizable way. Thus for N 
deltas we have 



E 

02x2 



2x2 



Ai 02x2 02x2 

02x2 



E 

02x2 



A2 02x2 

E A3 2X 2 
02x2 



... 



whith E and Aj defined in 



U 02x2 02x2 02x2 

and 



2x2 



02x2 

E A N 



V 



0. 



(21) 



2(AT+l)x2(AT+l) 



u 



1 1 





V 





-ika 



(22) 



We found the condition to be factorizable in terms of the functions hj(e) in a similar 
manner to that of the periodic chain. The eigenenergies of the system are the roots of: 



sin 



» ■ A(e;ax, ... ,a N ) 







(23) 



where 
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i iV even 

A{e; ai,...,a N ) = 2 N (h 1 • . . . • h N ) - 2 N ' 2 ^...(N - 2)...h k + 

p 

+ 2 N -*J^h i ...(N-4)...h k - ...(-l)f- 1 2 2 ^'^...(2).../ ifc + (-l)f (24) 



iV odd 



A(e;a h ... ,a N ) = 2 N ~\h l ■ . . . ■ h N ) - 2 N ~ 3 ^ ^...(N - 2)...h k + 

p 

+ 2 N - 5 J^h t ...(N-A)...h k - ...(.-l)— 2 2 ^'//,...C5)...//, + (25) 

P P 

TV— 1 

(-1)— {hi + h 3 + h 5 + ... + h N - 2 + h N ) 



and here the symbol hi...(M)...hk means the sum o/ a// possible products of M 

different hi 's with the following rule for each product: the first index has to be odd, the 
indices must follow an increasing order and to an odd index must always follow an even 
index and reciprocally. 

From (|23|) we see that e = nil ; n G N, are always eigenvalues of any finite length 
disordered chain wathever the species in it. 

2.3 Density of electronic states for an infinite chain 

The basic assumption concerning the electronic energies of a random chain is that the 
distribution of levels converges to a limiting distribution as the number of atoms goes 
to infinity, which is the same for almost all atomic sequences (except for a fraction that 
goes to zero as iV — >• oo) as long as the concentrations of the different species remain 
fixed. The existence of this property is needed by the thermodynamic limit: all quantities 
charaterizing macroscopically a random infinite chain with fixed species and concentra- 
tions cannot depend on the order of any finite-length piece of it. The method used to 
obtain the density of states in random chains is essentially due to James and Ginzbarg 
II and Schmidt 0] , who derived a functional equation which supplies the limiting distri- 
bution for the positive part of the spectrum. We have reconstructed and extended the 
method to provide the density of states for all energies ( Appendix O ). The functional 
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equation is, 

m 

W(y>) = $> {W (T-\ V )) - W (T^O)) } (26a) 

i=l 

W(y? + rvr) = W(p) + r (26b) 

W(?r) = 1 (26c) 

W(y?) is monotonically increasing in <£> (26d) 

where y> G (0, 7r] and m is the number of species composing the chain , jOj the ith species 
concentration and T~ i^p) are the functions, 

Tr\ V ) = arctan ^(e) - ^tt) (27) 

with hi(e) defined in (|). Solving (^) for different values of the energy, the density of 
electronic states per atom can be obtained from, 

c(e) = + (28) 

where K = -£™ lPi [W (3^(0))]. 

We solve this equation numerically. The range of <£> is discretized in p equally spaced 
points and W is represented by its values at those points which using a selfconsistent 
algorithm we have been able to calculate up to an error of 10~ 15 . Some examples of 
distributions of random chains are shown in Fig. 

The most interesting feature of the different spectra is the localized peak structure that 
appears for certain energy intervals. Those irregularities had been reported several years 
a §° [010 an d now the extremely accurate numerical algorithms dramatically confirm. 
Although the general interpretation of the behaviour of the distributions for these systems 
is non trivial one can say as a general rule (and following the suggestion of Agacy and 
Borland ||) that the peaked regions are remarkable in ranges of the spectrum which are 
forbidden for some of the species involved but not for all of them since in this latter case 
the range is also not allowed for the random system. A clear view of this explanation can 
be obtained by looking at a random chain of two species A and B. In a region forbidden 
for the A-type chain the allowed energies appearing in the spectrum are due to atomic 
clusters involving B atoms. As a result the A atoms surrounding these groups isolate 
those energies causing the density of states to decrease when we move in a tiny region 
around each of the levels and therefore giving rise to a fluctuating distribution. Thus 
one could reproduce the energies where the density of states would be more prominent, 
from the eigenvalues of certain atomic clusters in which the B species have a substantial 
contribution as shown in Fig. [|. In the (a) example it is clear how the eigenstates clusterize 
around the more peaked regions of the distribution. However one needs to consider a huge 
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(a) 2 species. (f-){pi} ■ -2{0.5}, -0.25{0.5}. The dotted line is the density of states for 
the ordered diatomic chain 




(b) 3 species. (%){Pi} : 10{0.3}, -1{0.3}, -2{0.4} 




(c) 8 species. (fJiPi} ■ -0.5{0.1}, -1{0.1}, -1.5{0.1}, -2{0.1}, -2.5{0.1}, -3{0.1}, -3.5{0.1}, 4{0.3} 



Figure 4: Distributions of electronic states for different configurations of the random 
chain. The thick line over the abscisa axis marks the common forbidden bands for the 
species composing the array. 5001 points have been taken for ip and Ae = 4 • 10~ 3 . 
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(a) pa= pb = 0.5 



B 

BB 

BBAB 
BAAB 



0.5 



1 1.5 

E 



2.5 



(b) PA = 0.9, p B = 0.1 



Figure 5: Density of states for a 2 species random chain with (^-J = —2 and (^-) = —0.25 
in a range forbidden for the A-type chain. The color vertical lines indicate the position of 
the eigenstates of different atomic B-groups surrounded by a certain number of A atoms. 
5001 points have been taken for (p and Ae = 2 • 10~ 3 . 
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number of clusters in order to reproduce all the maxima due to the equal concentrations 
of the species. Decreasing the B concentration we see in (b) how the sharp points can be 
quite easily predicted. 

Another feature that can be seen in the different representations in Fig. f|, is how 
the distribution approaches a smoother curve that resembles the one corresponding to a 
periodic chain as the energy grows: the greater the electron's energy the less it feels the 
presence of disorder. 

In Fig. |6], the evolution of the density of states can be observed for a pure one species 
chain as we dope substitutional^ with atoms of a different kind. 



2.4 The degree of localization of the electronic states 

Since Anderson || the disorder in the structures has been accepted as the main effect 
responsible of the spatial localization of the eigenstates of the systems. Just few years 
ago it was shown how if the disorder is correlated the localization can be suppressed for 
certain energies [@[10|. In 1999 the first experimental evidence that correlations inhibit 
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localization of states in disordered low- dimensional systems was reported fTTfl . 

Let us check the degree of localization for the electronic states in random chains with 
positive energy. What we specifically calculate is the average logarithmic decay per atom 
of the square of the envelope of the real wave function which seems an adequate quantity 
for such a characterization (12| . 

In the jth sector we write: 

= Aj cos (kx - k(j - l)a + fa) (29) 



and we look for: 



j=l \ J ' 

The connection equations of the wave function at the borders of the different sectors 
imply: 

4gi = 1 + ^ A^J cos 2 (A;a + fa) - (J^J sin(2fca + 20,). (31) 

Again, we will be using the functional equation to carry out the calculation so using the 
phase ipj as in the Appendix ||, it is just a matter of algebra to obtain: 

A) +l f 2 \ 2 sin 2 (£;a) 

m,aj) = —fir - + \J^.J 1 - 2 cos(fca) tan(^) + tan 2 (^) + 

/2 \ 2sin(fca)(cos(fca) -tan(y?j)) 2 
\ka,j J 1 — 2 cos(&a) tan(yjj) + tan 2 (c,Oj) 



and the same arguments apearing on the Appendix lead us to write ( p0|) as 

(log JO = f> r^^-logF^a^dcp. (33) 

Integrating by parts and using the equations for W(y>), we finally obtain: 

(logF) = f^jftlogFfoo,) - f> fw(y) 1 ^J'^ dy. (34) 

In the Fig. [7] some examples of the degree of localization as a function of the energy 
are shown. As can be seen the degree of localization has a tendency to decrease in the 
peaks of the spectrum and increase in the troughs, and for ka — rm ; n G N the states 
are always extended for any chain. This last result had been already obtained using the 
trasfer matrix technique ITS for finite chains. So the system with completely uncorrelated 
disorder can support extented states. In fact an infinite number of isolated resonances is 
present altough mobility edge for the electrons does not exist and therefore one cannot 
speak of an Anderson transition in a strict sense for this particularly simple model. 




Figure 7: Degree of localization for the electronic states (black line) and density of states 
(grey line) for a random chain. 5001 points have benn taken for ip to solve the functional 
equation and A(ka) = 4 • 10" 3 . 
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3 Fractality 

Let us observe the peaked regions of the distributions. At first sight the irregularities 
make the density of states appear hardly different iable within those intervals. Would 
this pattern still hold as we look deeper into the distribution? Is the distribution non 
differentiable? In order to answer these questions one has to improve the numerical 
calculation to be able to see the real density of states in smaller energy ranges. In Fig. |8] 
the spectrum of a certain random chain is shown for shorter and shorter energy intervals. 
As can be seen as the energy domain is made smaller and thus increasing the accuracy 
of the numerical algorithm the distribution reveals a finer structure: new sharp points 
appear and the density does not evolve smoothly. We have been extremely careful with 
the numerical algorithms and we are pretty sure that the observed irregularities are not 
due to numerical errors. For the representations in Fig. |^ we have proceeded doubling the 
number of points taken for (p in order to represent the functional equation and checking 
the convergence of the density of states at each step until the desired accuracy ( in all 
cases the average variation of G(e) in the last step relative to its domain was less than 
0.75% ). The final parameters were: 



points for (p Ae 



Fig. 


|(a) 


5001 


7.5- 


10" 4 


Fig. 


|(b) 


35001 


2.5- 


10~ 5 


Fig. 


|(c) 


150001 


2.5 • 


■ 10~ 6 



As a final check we repeated the procedure for the same random chain, with the same 
paremeters but in an apparently smooth region of the density of states obtaining the 
results shown in Fig. |{J From this we can see that the irregularities in certain intervals of 
the spectrum are not due to errors of the numerical computation. 

The numerical calculation of the density of states suggests the possibility for the distri- 
butions to show a quasi-fractal behaviour in some energy ranges: non differentiability for 
any point inside these regions and irregular aspect whatever the scale. The differentiabil- 
ity of the distribution is a consequence of a regular (almost homogeneous) distribution of 
the states inside a small energy interval. A clear view of the differences between the way 
in which states appear inside an irregular region and a smooth region of the density can 
be obtained representing the energy spaces for adjacent levels. In Fig. [Tl] these spacing 
distributions are shown for the first 100 levels that appear inside an irregular zone and a 
smooth one. Increasing the length of the random chain we see how the spacing distribu- 
tion for the levels in the smooth region becomes more and more homogeneous as the first 
100 levels are included in a smaller energy interval while the spacings for the levels in the 
irregular zone does not show a defined tendency nor a homogeneus distribution. In fact 
these last spacing distributions exhibit the same aspect whatever the scale would be. 
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Figure 9: Density of electronic states in different energy ranges for a two species random 
chain (^-) = — 2 and (^) = —0.25 with equal concentrations smooth region). 
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The fractality of the density of states does not depend on any particular parameter of 
the random chain; its presence is a consequence of the disorder and the dimensionality of 
the system and it is in this sense a universal effect. This fractal behaviour might be related 
to the fractal conductance fluctuations observed in gold nanowires 0] in the mesoscopic 
regime suggesting a connection between disorder and coherent transport. 



1200 deltas 




4800 deltas 



9600 deltas 






Figure 10: Energy spacings for the first 100 eigenvalues appearing inside an irregular zone 
( e > 1.4, red ) and inside a regular one ( e > n, blue ). The point n on the abscisa axis 
represents the energy distance between the (n — 1) and n level. All the sequences are 
composed of the species (^-) = —2 and (^-) = —0.25 with equal concentrations. For a 
fixed length all the states calculated correspond to the same random sequence. The states 
have been calculated by finding the transmission of the wave function through the chain 
with fixed end-points boundary conditions. 



4 Concluding remarks 

To summarize we would like to point out the main results we have obtained in this paper. 
First, we emphasize that the band structure provided by (|15|)-(|TTD is not just exact but 
also extremely useful from the point of view of computer algebra calculations. In fact we 
have carried out various profiles for the curves provided for these conditions until N=30 or 
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more using just few seconds of a lap-top regular computer. The reason for that lies mainly 
in the systematic use of the form, products and combinations of the /i(e)-function defined 
by The Luttinger theorem has been generalized to an arbitrary number of different 
delta potentials and the Saxon-Hutner Conjecture has also been definitively established. 
Similar expressions have been obtained for the case of a disordered finite chain and the 
extended functional equation method has been used to look specifically to the density 
of electronic states for a random infinite array. By obtaining the degree of localization 
of the electronic states in the random chains, we have a more complete picture of the 
role the levels play in the transport processes and the existence of extended states under 
uncorrelated disorder has been confirmed. The results are not only in agreement with what 
was expected in more sophisticated models but also go beyond them and we can account 
for universal properties of the transport effects by detecting fractality in the conductance. 
This work is far from being finished as we want to ascertain ourselves whether this effect 
can be measurable in terms of the fractal dimension and critical exponents. At the 
moment, however, we can already offer a large bundle of properties that arise from a 
model whose simplicity is not only a shortcoming but rather an advantange to study 
complex properties in a benchproof easily manageable. 



A Appendix: Gaps theorem 

Let us consider a finite chain of N delta potentials of different species with fixed end- 
points boundary conditions ( Fig. [TTJ ). Inside the jth sector we write the wave function 
of the electron 

^j(x) = A,,' 1 ''' + B je - ikx \ (35) 

where Xj = x — (j — l)a. 



V(x) 



N+l 



x=(N+l)a 



Figure 11: Deltas chain 
Imposing the conditions at the border x = ja one can obtain the relationship between 
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the amplitudes of adjacent sectors: 

i ^ ^ika —i ^—ika ^ / A \ 

(36) 



A j+1 \ {±-ity ka ^- %ka \(A i 



M 

and M is the transmission matrix of the sectors. The boundary conditions are: 

A x + B x = (37a) 
A N+l e ika + B N+1 e~ lka = 0. (37b) 

A ■ 

Let us define inside each sector Tj = —ie tka —^-. The Eqs. ( |37|) become: 

n = ie ika (38a) 
r N+1 = ie~ ika , (38b) 

and from ( |36| ) the transmission for r,- is 

CjTj + 1 



^ = f^r (39) 
'i c i 



where = (i — ka j )e lka . Let us use z,- = The boundary conditions now read: 

3 y 31 3 l+iTie lka 



z x = (40a) 
Zn+i = oo- (40b) 



A straightforward calculation using (|3~9"D lead us to: 



= o, : ■ ( 41 ) 



2hj(ka 



z 



with hj(ka) = cos(ka) + [Jjzj — ^ • Now let us consider the values of ka satisfying 

\hj(ka)\ > 1 for every aj in the chain (i.e. the common forbidden bands for all the one 
species chains). As \z\\ < 1 it is easy to proof by induction using (f4~T|) that \zj\ < 1 

< 1. So we have |£jv+i| < 1 for arbitrary N and the final boundary condition cannot 
be satisfied whatever the length of the chain. Therefore those ka values are not allowed in 
the mixed chain. Notice that the order in which the different species appear is irrelevant. 
The same conclusion holds for both periodic and disordered chains. 

The result is also valid for complex values of k (i.e. negative energies). 
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B Appendix: Functional equation 

The limiting distribution of electronic states for a random chain composed of m different 
species each one with concentration pi, will be given by the average of the spectra for all 
sequences which have the given concentrations in the limit as the number of atoms goes 
to infinity. 

Let us consider the same chain as in the Appendix [A| and the quantities Zj. Since 
z\ = and the transmission (^TJ) is real every Zj will also be real and we can define a 
phase in each sector tpj so that tan(^) = Zj, which yields the result: 



<p j+1 = Tjfa) = arctan - ^ j . (42) 

We need the transmission of the phase to be an increasing continuous function. Thus we 
define: 

f arctan ( , , > \ — -, — r ) si ip« G ( — arctan (2hj(e)) 

\2 hj (e)-t W {n)J \ K A ' J \ (43) 

arctan ( OL — s ) + rr si (fj G ( arctan (2hj(e)) , £ 



^•(^ + rwr) = 7}(^) + rar ^ G | . (44) 

This can be easily written by using the inverse function as: 

7- te ) = arctan (2^) - j-j-j) (45) 

T^ 1 ^- + n?r) = + n?r tpj G (0, tt] . (46) 

The boundary conditions are: 

<Pi = (47) 

7T 

Pn+i = — + nn n G Z. (48) 

Once the transmission of the phases can be written in a uniquely way, 93^+1 (e) is a 
continuous function of e. Let us assume the latter to be an increasing function. Then 

if €\ < €2 the quantity (| ° Af + 1 ^ 2 ) ( Pn+i(^i) . g ^ e number of times the final boundary 

71 

condition has been satisfied from ex to ti with an error smaller than 1 and therefore it 
also represents the number of eigenstates in the interval t\ < e < e 2 . So it is clear that 
we can write the density of states per atom of the chain whatever the behaviour of the 
function ip N+ i(e) would be as: 



1 



<PN+l{t+de) _ i/?iy+i(e) 
N N 



G < £ > = H * (49 > 
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As ( f47| ) holds for all e, we can write 



yjv+i(e) . 



i=i j i 



(50) 



for each value of e. Thus y>iY t r 1 ^' 1 is the average over all atoms of the advanced phase which 



A' 



we denote (Aip)(e). As Tj((fj) — <pj is a periodic function with period it the average (Aip) 
can be calculated using the distribution functions of the phases <pj( mod n). 



(51) 



where Wj(<£>) with <£> G (0, 7r] is the probability that (pj( mod 7r) lies in the interval (0, ip\. 
Also we impose: 



Wj(ip + m) = Wj(ip) +r ip E (0, 7r] . 



(52) 



The probability distribution for ipj depends on the distribution for (pj-i: ipj(mod tt) 
lies in (0,<^] if and only if <fij-i( mod tt) lies in (7^(0), T~ x (</?)] ( mod 7r). There exist 
integers s and t such that 



^- 1 (0)(mod tt) = Tj\0) = Tr\0) - S n 
Tf 1 M( mod tt) = TT 1 ^) = t;- 1 ^) -«7T 



and < Tj~ (ip) — 7^ _1 (0) < 7r from ( f45"|) thus £ must be equal to s or s + 1. 
1. |t=s 



(53) 
(54) 



In this case (Fig. 12): 

Wfa) = W 3 . x (TJ 1 ^)) - W^_i (T7 X (0)). 



Tj(0) Tj(9) 

' i 



t, to) t; (<p) 



4- 



Figure 12: Example t = s = 2 



2. |t=s+l 

In this case (Fig. [TT3| ) : 
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Tj'(<p) T(o) 






















71 . ■ ' • . v * ' 


1 

271. " 




3Jt.-*' 



Figure 13: Example s = 2 and £ = 3 



Using (HD we can unify both cases as = Wj-i (Tf 1 ^)) - V^_i (7} _1 (0)). Thus 

the distribution functions are the solutions of the equations 

Wjfr) = Wj-i (Tr\ip)) - Wj-i {Tf\0)) (55a) 
Wj((p + rvr) = Wj(ip) + r (55b) 
Wjiir) = 1 (55c) 
Wj(ip) is an increasing function of ip (55d) 

with ip G (0, 7r]. 

The random chain is composed of m different species with lengths a%, . . . , a m in con- 
centrations Pi, . . . ,p m (Y2iLoPi = -*-)■ We nave now to take into account all possible 
sequences which have the given concentrations. Thus we only know the probability of 
finding certain species in the position x = ja and therefore we must take the average. 
The Eq. (^l|) becomes: 

and the equations for the distribution functions become: 

m 

Wjitp) = $> {Wj-i (TTH<P)) - Wj-! {Tr\0))} (57a) 
i=i 

Wj{<p + rw) = Wj(<p) + r (57b) 

Wjin) = 1 (57c) 

Wj(tp) is an increasing function of <p (57d) 

with p G (0, 7r]. Notice that now is the probability that pj( mod 7r) lies in the 

interval (0, <p] in some of the possible sequences. Let us write ([56]) as 

(A^ = E^/ — r^ 1 ~ <P} *p (58) 
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1 N 

where W N {if) = — Wj((p), which is the average distribution function of the chain. 

3=1 

Let us take the limit N — > oo and denote W(t^) = lim W N (ip). It is not hard to see 

N—>oo 

that this function will satisfy: 

m 

Wfo>) = ^ftlWft" 1 ^)) -W(7T- 1 (0))} (59a) 
i=i 

W(y? + rvr) = W(<p) + r (59b) 
W(?r) = 1 (59c) 
W(y?) is an increasing function of (59d) 

with (y9 G (0,7r]. And, 

( = e ft r {^(v) - *>> ^- (so) 

Schmidt has proved 0] that the solution of Eqs. fl59"D for each energy is unique and 
continuous. To carry out the integration in (|60|) we use the existence of a value <po — | 
such that 7i(<^o) verifies Ti(ip ) = ipi = n for any species c^. 

= Eft / — % - = Eft / — % -v}^, 

t! Jo dtp ^ J VQ dip 

integrating by parts 



(A?) =Eft w (v) - ^} - Eft / wM-p# + 

i=l Po i=l ■'Vo a( P 

+ Eft / w (^)^ 



=1 • / ¥>0 

and using ( |59a| ) 

W ^"£^ = 5> / 



Eft / W(^)^p^ = Eft / w (T^)) ^ = 

r +W W(9)d6 + 7rf>W (T'-^O)) . 

>^ 031 • i 
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Going back to the expression for (Aip) we obtain: 

W{ip)dip- \ 

i—l •>•£() Jfl 

From ( |59b|) is easy to see that 



/ W((p)dtp- I W(9)d8 = ipo-^ 



and finally the following relationship holds for a certain value e of the energy: 



(Ar:) -/C = -f; ft W(7-- 1 (0)). (61) 



7T 

i=l 



Solving the Eqs. for different values of e one can calculate the density of states of 
the random chain from 

G(£) = \^±f)_^M. (62) 
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